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^ ' ABSTRACT 

, ^ ' Context. The precise and accurate modelling of a terrestrial planet like Venus is an exciting and challenging topic, all the 
more interesting since it can be compared with that of the Earth for which such a modelling has already been achieved at the 
milliarcsecond level 

Aims. We want to complete a previous study (Cottereau and Souchay, 2009), by determining at the milliarcsecond level the 
polhody, i.e. the torque-free motion of the axis of angular momentum of a rigid Venus in a body-fixed frame, as well as the 
nutation of its third axis of figure in space, which is fundamental from an observational point of view. 

Methods. We use the same theoretical framework as Kinoshita (1977) did to determine the precession- nutation motion of a rigid 
• , Earth. It is based on a representation of the rotation of rigid Venus, with the help of Andoyer variables and a set of canonical 
' equations in Hamiltonian formalism 
Results. In a first part we have computed the polhody, i.e. the respective free rotational motion of the axis of angular momentum 
^ . of Venus with respect to a body-fixed frame. We have shown that this motion is highly elliptical, with a very long period of 525 
' cy to be compared with 430 d for the Earth. This is due to the very small dynamical flattening of Venus in comparison with our 
C/3 planet. In a second part we have computed precisely the Oppolzer terms which allow to represent the motion in space of the 
■ third Venus figure axis with respect to Venus angular momentum axis, under the infiuence of the solar gravitational torque. We 
' have determined the corresponding tables of coefficients of nutation of the third figure axis both in longitude and in obliquity 
, due to the Sun, which are of the same order of amplitude as for the Earth. We have shown that the coefficients of nutation 
^ ■ for the third figure axis are significantly different from those of the angular momentum axis on the contrary of the Earth. Our 
\^ analytical results have been validated by a numerical integration which revealed the indirect planetary effects. 
CNj , Conclusions. This paper is a complementary study of Cottereau and Souchay (2009). It gives a precise determination both of 
' the torque free motion of Venus, as well as the nutation of the third Venus figure axis in space for a short time scale, when 
' considering the planet as a rigid body. 

^ ; 

' Key words. Venus rotation, rmtation 
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J> 1. Introduction We have performed a full calculation of the coefficients 

of nutation of Venus and presented the complete tables of 

Venus which can be considered as the twin sister of the nutation in longitude A* and obliquity Ae, for the axis of 

\ Earth, in view of its global characteristics (size, mass, den- angular momentum due to both the dynamical flattening 

sity), has been the subject of a good amount of investi- and triaxiality of the planet. 

gations on very long time scales, to understand its slow In this paper, the study begun in Cottereau and 

retrograde rotation (243 d) and its rather small obliquity Souchay (2009) is completed. First in section [5] we con- 

(2°. 63) (Goldstein 1964; Carpenter 1964; Goldreich and sider the torque free rotational motion of a rigid Venus. 

Peale 1970; Lago and Cazenave 1979; Dobrovoskis 1980; We recall the parametrization of Kinoshita (1977) and the 

Yoder 1995; Correia and Laskar 2001, 2003). Habibullin important equations of Kinoshita (1972) which are used 

(1995) made an analytical study on the rotation of a rigid to solve this torque free motion. The important charac- 

Venus. In Cottereau and Souchay (2009) we presented an teristics (amplitude, period, trajectory) of the free motion 

alternative study, from a theoretical framework already are given . Cottereau and Souchay (2009) supposed that 

used by Kinoshita (1977) for the rigid Earth. We made an the relative angular distances between the three poles (of 

accurate description of the motion of rotation of Venus at angular momentum, figure and rotation) are very small as 

short time scale. We calculated the ecliptic coordinates of it is the case for the Earth. In this paper we want to deter- 

Venus orbital pole and the reference point jov which is the mine accurately the motion of the third Venus figure axis 

equivalent of the vernal equinox for Venus. Our value for which is the fundamental one in an observational point of 

the precession in longitude was ^ = 4474" .35t/cy±66.5. view. To do that we reject the hypothesis of coincidence 
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of the poles. Thus we determine the Oppolzer terms de- 
pending on the dynamical flattening and the trixiality of 
Venus. Then we compare these terms with the correspond- 
ing terms of nutation for the axis of angular momentum, 
as determined by Cottereau and Souchay (2009) and the 
Oppolzer terms determined by Kinoshita (1977) for the 
Earth in section [31 We give the complete tables of the 
coefficients of nutation of the third figure axis of Venus. 
We compare them, with the coefficients of nutation of the 
angular momentum axis, taken in Cottereau and Souchay 
(2009) (section 0]) . Finally in section [5] we determine the 
nutation of the angular momentum axis by numerical inte- 
gration using the ephemeris DE405. We validate the ana- 
lytical results of Cottereau and Souchay (2009) down to a 
precision of the order of a relative 10~^. We show that the 
discrepancies between the numerical integration and ana- 
lytical results (Cottereau and Souchay, 2009) are caused 
by the indirect planetary effect, i.e the small contribution 
to the nutation, which is due to the periodic oscillations 
of the orbital motion of Venus. In this paper as it was 
the case in Cottereau and Souchay (2009), our domain of 
validity is roughly 3000 years. 



2. Torque free motion for rigid Venus 

2.1. Equation of torque-free motion 

We consider the problem of the rotational motion of the 
rigid Venus in absence of any external force. We note (0, 
X, Y, Z) the inertial frame and (0,x,y,z) the cartesian co- 
ordinates fixed to the rigid body of the planet (see FiglT]). 
The orientation of Venus, with respect to the inertial axes, 
is determined through the Euler angles : hf, If , cj). The pa- 
rameter (j) gives the position of the prime meridian (0,x) 
with respect to 7oy (Kinoshita, 1977). The angular mo- 
mentum axis (hereafter denoted AMA) of Venus is the 
axis directed along G. The components of the angular mo- 



plan normal to G 




inertial plane of Venus 
at epoch J2000.0 



equatorial plane to Venus 



Fig. 1. Relation between the Euler angles and the 
Andoyer variables. 



mentum referred to the body-fixed axes are : 

Lx = AlUx = A(hf sin If sin (p + If cos (f) 
Ly = Buy ~ B{hf sin// cos<j} — If saKp) 
Lz — Cujz = C{hf cos If + (p) 



(1) 



where A, C are the principal moments of inertia of 
Venus. Moreover the kinetic energy is : 



(2) 



To describe the torque free rotational motion we use the 
Andoyer variables (Andoyer, 1923; Kinoshita, 1972) (see 
FigUl): 

— L the component of the angular momentum along the 
Oz axis 

— H the component of the angular momentum along the 
OZ axis 

— G the amplitude of the angular momentum of Venus 

— / the angle between the origin meridian Ox and the 
node P 

— h the longitude of the node of the AMA with respect 
to 7oy 

— g the longitude of the node of the plane (0, X, Y) with 
respect to Q and to the equatorial plane. 

This parametrization is described in details in Kinoshita 
(1972, 1977) and in Cottereau and Souchay (2009). From 
these definitions we have : 



L — G cos J, H = G cos / 



(3) 



where /, J are respectively the angle between the AMA 
and the inertial axis (O, Z), and the angle between the 
AMA and the third figure axis (hereafter denoted TEA). 
Using spherical trigonometry, we determine the following 
relation between the variables: 



(4) 



The components of the angular momentum vector with 
the Andoyer variables are : 



Lx 



\/G'2 - 7.2 sin/ 
L. 



L2, 



(5) 



The Hamiltonien for the torque-free motion of Venus cor- 
responding to the kinetic energy is : 



2^ A 



(6) 



The Hamiltonian H does not depend on the time and 
it is free from g and h. Thus the number of degrees of 
freedom of the torque free motion is one. Deprit (1967) 
has characterized this motion by studying the isoenergetic 
curves in the phase plane L-1. However we prefer to study 
the isoenergetic curves in the (J, I) phase space using the 
equation ([3]). This allows a good description of the position 
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Circulation motion 



Libratlon motion 



Fig. 2. Isoenergetic curves of Venus in the (J, /) phase 
space of the torque free motion. Two motion are showed: 
" the hbration motion" and the " circulation motion" . 

of the AM A with respect to the TFA. Fig H] shows the two 
possible motions : one is the motion of libration and the 
other one the motion of circulation. Using the values of the 
relative differences of the moments of inertia of Williams 
(private communication): 

^ = 5.519x10-^ = 3.290 x 10- 

B — A 

= 2.228 X 10— 

Mi?2 



(7) 



and addopting 



C 



c 
WW 



1.000016. 



0.3360 (Yoder, 1995), we get 
= 1.000010. We introduce here j which 
corresponds to the minimum value of J when ^ = f ■ We 
denote by ji the value reached by j on the separatrix. The 
libration motion is not possible if J < ji or J > (n — ji). 
The numerical values of ji can be determined using the 
equations of Kinoshita (1972). We have : 



ji = arccos 
with 

and 

1 _ 1 



2e 



1(1 



A' 



(8) 



(9) 



(10) 



where e measures the triaxiality of the rigid body 
(Andoyer, 1923). Applying these formulas to Venus we 
get : ji = 52°. 23. 

2.2. Canonical transformations 

To determine the torque free motion of Venus, we use the 
method described in full details by Kinoshita (1972). We 
only give here the important results and equations which 



are needed to apply this study to Venus. The Hamiltonian 
for the torque free motion is given by To solve the 
equations of motion, we perform a canonical transforma- 
tion which replaces the angular momentum by its action 
variable (Goldstein, 1964). First we make two intermediate 
transformations to simplify the computation. 

2.2.1. First transformation 

We change (L, G, 1, g) to (ai, 0:2, wi, U2) using 
the Hamilton Jacobi method (Chazy, 1953). We note 
S{L,G,ai,a2) the characteristic function. We have the 
Hamilton- Jacobi equation : 

_ 1 sin^Z cos2/ 95 ,dS^ 1 55 2 .... 



where 



ai ^ -to, ^2 = = 

09 




7^ — eal cos 21 

— —dl + a2g 

1 — ecos2/ 



(12) 



(13) 



(14) 



The constants e and D are given by © and ([TU]) . Because 
the new Hamiltonian depends on only one of the momenta 
ai, we have : 



Ul 



t + fii 



dS _ dS_D_ 

dai 97 7 



U2 = P2 



dS 
da2 



(15) 



2.2.2. Second transformation 



Then we make another transformation which changes 
(ai,a2, wi,U2) to {L,G.,l,g). The two momenta are de- 
fined as follows : 



L = 7 

G — a2 — G 



(16) 



With these new coordinates the Hamiltonian Fq becomes 

Fo = l^LV + i(i + >V (17) 

Thanks to the Hamilton equations the conjugate variables 
are: 



Using (ITS]) we get J from: 

L — G cos J 



(18) 
(19) 

(20) 
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and 



cos J = yi — {1 + e) sin j — cos j y 1 — e tan'^ j (21) 

where j is defined in the subsection 12.11 Now we can in- 
troduce the canonical transformation which replaces the 
angular momentum by its action variable. 



2.2.3. Action variables 

The action variables are given by: 



2.3. Venus free rotation 

As a result we can use a development of our variables 
g, I and J with respect to j and e. The results of these 
technical developments (Kinoshita, 1972) are summarized 
in the appendix. As it is the case for the Earth we can 
suppose that the angle j of Venus is very small. Table [T] 

Table 1. Important values of the free motion of Venus 




L 



G = G = G = a2 



eG2 cos 21 



Constant numerical values 



e cos 21 



dl 



(22) 

/ and g are canonically conjugate variables respectively 
to L and G. We make the following transformation to 
simplify the calculation: 



cos 25 



cos 21 



or cos 21 



cos 2(5 



1 — ecos2Z l + ecos2(5 

We obtain : 

G r - k'^cos^s 



(23) 



c 

6 

B 


1.000016 




1.000010 




e 


0.230769 




c 

D 


-0.000013 






0.0119495 


rd/cy 


Ti 


525.81 cy 




Ug 


0.0258549 


rd/d 


Tg 


-243.02 d 





^ 2. 
= GAo, 
where 

and 

Ao = 
with 



V6 



1 -1- e cos 25 



dS 



1 



(6-1) &=(cos2j)- 



(24) 



(25) 



gives the numerical values of the important constants used 
in the theory of the free motion of Venus. To calculate 
these values we take jf^ (Yoder, 1995) and the values of 
the moment of inertia given in ([T]) . The value of D and b 
are unknown. Finally we find: 



1 = 1* - 0.0576923j" 
G 



sin 21 



1 + e cos 21 



0(/) 



E{k)F{x, k') + K{k)E{x, k') - K{k)F{x, fc')]26) ^ = 5 + [ - 0.973009(r - f) 



X 



sm 



k' = \J\ - fc2. 



-^0.307692 X j2 [0.230769 



sin 21 



(27) 



1 -I- e cos 21 



In equation (051) Kik) is the complete elliptic integral of 
the first kind with modulus fc^, i?(fc) that of the second 
kind. i^(x, fc') is an incomplete elliptic integral of the first 
kind, -E(x, fc') that of the second kind, Aq is a Heumann 
lambda function (Byrd and Friedman, 1954). The com- 
plete solution of being of no interest here, we only 
give the mean motion of the angular variables / and g that 
will be used hereafter (Kinoshita, 1972, 1992). Thanks to 
the canonical transformations the time variation of / and 
g can be determined. We have : 



J 



2.05548(/* - 0]] 
= j\l 1.28713 + 0.287129 cos 2[-f 0(e^ j^) 



g = Jigt + 132 = ngt + go 



I = nit 
where : 

ni 



ttG 
2KD 



cos J 



G 



-Gi^-j)sm'j. 



(28) 
(29) 

(30) 
(31) 



tan/* = 0.797007 tan/. 
Using the developments for a small e, we get: 

/ = /~- 0.111538 sin2/~-|- 0.006221 sin4/ -I- 0(e^j2) 

g = g + 0.111538 sin2/- 0.006221 sin4/ -f 0(e^ j^) 
J= J[l. 00933 + 0.11154 cos2/] -I- 0(e^/) 

J = j + 0.115385 tanj 
+0.05325 tan j(- + — tan^ j) + O(e^). 



(32) 



Here the epoch of time t is defined such that / = ^ 



at 



t = Q. Now we can give the development of the variables 
g and / which describe the torque free motion. 



The main limitation of our calculation being the uncer- 
tainty on the ratio -j^^ of Venus (Yoder, 1995), our poly- 
nomial expansions must be done accordingly. 4*'' order 
terms are too small compared to our level of accuracy and 
have been discarded. Projecting the pole of Venus on the 
(X, Y) plane with a value of j = 0.01 rd for the develop- 
ments above, we can plot the free motion of Venus (Fig 1^1). 

We note here that the value of j has been chosen ar- 
bitrarily, for it does not affect in a significant manner the 
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Fig. 3. (X, Y) free motion of Venus in space for five hun- 
dred century time space (red and blue curves). The green 
curve represents a circle with radius of r — j. We take 
j = 0.01 rd. 

polhody, except the amplitude. We see that the torque 
free motion of Venus is an elliptic motion, as it is the case 
for the Earth. The rotational free motion of Venus, with a 
period T/ = 525 centuries (Tab[T]) is much slower than the 
Earth one (303 d). If we consider the elastic Earth (ocean 
and atmosphere), the torque free motion has a period of 
432 d, significantly larger than in the rigid case. The at- 
mosphere of Venus is much denser than that of the Earth. 
So it would be interesting to study the torque free motion 
of the elastic Venus in a next paper. 

3. Rigid Venus forced rotational motion 

In Cottereau and Souchay (2009) we supposed that the 
relative angular distances between the three poles of Venus 
(pole of angular momentum, of figure and of rotation) are 
very small, as it is the case for the Earth. In this section 
we want to determine the motion of the TEA, which is the 
fundamental one from an observational point of view. To 
do so we reject the hypothesis of coincidence of the poles. 
Using the spherical trigonometry in the triangle (P, Q, R) 
(see FigllJ) we determine the relations between the TEA 
and the AMA. Supposing that the angle J between the 
AMA and the TEA is small we obtain: 

sing + 0(j2) (33) 



sin/ 



0(J2 



(34) 



If — I + J cosg 

where / characterizes the obliquity and h the motion of 
precession- nutation in longitude of the AMA of Venus, hf 
and // correspond to the same definitions as h and /, but 
for the TEA instead of the AMA. This yields (Kinoshita, 
1977) : 

Ahf = Ah + A{:^^)+0{J^) (35) 



AIj 



AI 



sin/ 

A{J cos g) + 0{J^) 



where Ah and AI represent the variation of the nutation 
of Venus, respectively in longitude and in obliquity. The 



second terms in the right hand side of (|55)) and are 
the so-called Oppolzer terms. They represent the differ- 
ence between the nutation of the TEA and the nutation 
of the AMA (in longitude and in obliquity). These terms 
represent the differential effects of solid body tides on both 
axes. Developing the equations and ([55)) we obtain: 

;{AJ sing + JAg cosg 



A(:^) = 



sin / sin / 
J sing A/ cos/ 



/ 



(37) 

JsmgAg. (38) 
are functions of Ag and AJ. 



sm 

A(Jcosg) ~ A J cosg 

We see that ^ and 
We remind here that g and J caracterize the motion of 
the AMA with respect to the TEA. As the Earth has a 
fast rotation the TEA and the AMA can be considered 
identical. Venus has a slow rotation, so it is interesting to 
see what difference will emerge on the motion of the TEA. 
To determine Ag and A J we must solve the equations of 
motion. 

3.1. Equations of motion 

The Hamiltonian related to the rotational motion of Venus 
is (Cottereau and Souchay, 2009): 

K'' = Fo + E + E' + U (39) 

where Fq is the Hamiltonian for the free motion, E + E' is 
a component related to the motion of the orbit of Venus, 
which is caused by planetary perturbations. The expres- 
sion of Fo has been set in the previous section. The ex- 
pression oi E + E' is given in detail in Cottereau and 
Souchay (2009). U is the disturbing potential due to the 
external disturbing body considered. Here the sole exter- 
nal disturbing body is the Sun (The perturbation due to 
the planets can be neglected in first order), and its dis- 
turbing potential is given by : 



GM'2C~A 



_ _n A — ±j 

U = — P2{smd) + — ^P|(sin(5)cos2a] 

(40) 

where G is the gravitationnal constant, M' is the mass 
of the Sun, r is the distance between the Sun and Venus 
barycenters. a and 6 are respectively the planetocentric 
longitude and latitude of the Sun, with respect to the 
mean equator of Venus and with respect to a meridian 
of origin (therefore a must not be confused with the usual 
right ascension) . The P™ are the classical Legendre func- 
tions given by : 

(-1)™(1 - j;2)f (i"+™(a;2 - 1)« 



B 



(41) 



2"n! d"+™a; 
In (|40p we only consider the potential in first order. The 
method for solving the equation of motion is described in 
Kinoshita (1977). Only the final results are given. We have 



(36) Ag 



AJ 



1 , ^ .dWi 

1 1 dWi 
G^sinJ~dr 



cot/^) 



cot J 



dg 



(42) 
(43) 



6 

where 
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and 

1 



GM' 



A-B 



2C-A-B 



Pn sin S cos 2a 



dt. 



(44) 



We use a transformation described by Kinoshita(1977) 
and based on the Jacobi polynomials. It expresses a and 
S as functions of A and /3, respectively the longitude and 
the latitude of the Sun with respect to Venus mean orbital 
plane : 



P2(sin(5) =-(3cos2j-l) 



(Scos^/- l)P2(sm;9) 



--sin2/Pj(sin/3)cos2(A- h) 



sin 2 J 



— - sin 2/F2(sm /3) cos g 



- i(l + ecos/)(-l + 2ecos/) 

e=±l 

P2 {sin 13) sin(A — h — eg) 

— ^ if sin /(I + e cos/) 



e=±l 

32 



F2 (sin/3) cos(2A - 2h - eg) 



sin^ J 



- sin^ /P2(sm/3) cos 25 H — esinl 



e=±l 



(1 + e COS /)Pj (sin (3) sin(A - /i - 2eg) - 



16 



(1 + ecos/)2p|(sin/3)cos2(A -h- eg) 



e=±l 



(45) 



|(3cos2/- l)P2(sin/ 



P| (sin 5) cos 2a = 3 sin^ J 

cos2/ + i Y sm2I P^ {sin 13) sm{X- h-2el) 
1 



e=±l 



■ sin^ /P|(sin P) cos 2{X - h - el) 
o 

' ^ psin J(l + pcos J) 
p=±i 

3 

- - sin2/P2(sin/3) cos{2 pi + g) 

- Y ^(1 + ecos/)(-l + 2ecos/) 
£=±1 ^ 

P2 (sin /3) sin(A -h- 2pel - eg) 

— ^ -esin/(l + ecos/) 

e=±l 

P2^(sin 13) cos(2A - 2/i - 2peZ - eg) 



To simplify the calculations, we study separately the sym- 
metric part of (j44l) depending on the dynamical flattening 
and the antisymmetric one depending on the triaxiality of 
Venus. So the coefficient of the dynamical flattening will 
be noted with a "s" index and the coefficient depending 
on the triaxiality with an "a" index. From its definition 
above, we can set: /? ~ 0, for the latitude of the Sun with 
respect to Venus mean orbital plane can be considered as 
null. 

3.2. Oppolzer terms depending on the dynamical 
flattening 



Using the equations (|44)) and (|45)) . with /3 = 0, we have : 
GM' 2C-A-B 



Wsl = 



■ J[{^)'P2{sinS)]dt 



i(3cos2 J - 1)1^,10 - ^ sin2JW^,ii 



1 



■sin JWsi2 



(46) 



where 



dt 



— sin^/ / cos2{\- h){-)^dt] 



(47) 



Wsll 



X,.[-isin2/ j i(^)3cos(7 dt 
i(l-cos/) / cos(2A- 2/i-|-5)(-)^dt 



i(l + cos/) / cos{2\ - 2h - g){-)''^dt] 



(48) 
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Wsi2 = Ks' [- sin^ /i J{^f cos 2g dt 

-cos If j cos 2{X - h + g){^fdt 
-i(l + cos/)2 f cos2{X- h-g){^fdt] 



(49) 



and 



3GM' 2C - A~B 



(50) 



Thanks to the Hamihon equations, we obtain 

A Ir ^ 9 T cos2JcosJ^^^ 
As5 = —[-3 cos- J Wsw -T—i ^sii 



sin J 



1 



COS^ J Wsl2] 



cos I Ah + 0{J^ 



(51) 



We suppose that the angle J is smah as is the case for the 
Earth. This yields : 



A,5 = ^[-3Wso-jW,ii + ^Wsi2] 
-cos / A/i + 0(j2). 

We have also : 

AsJ = —[cos^ J -sin2J 

G og 8 



(52) 



dg 



1 dWsu 
G dg 



0{J). 



(53) 



Using ([5^ and we obtain the Oppolzer terms depend- 
ing on the dynamical flattening : 



, J sin g 



1 



sin/ Gsin/ 
-W,iiCos.g] + 0(J) 



■smg 



As (J cos g) = A( J) cos g — JAg sin g 



-cosg + Wsii sing] + 0{J). 



(54) 



(55) 



To solve the equations and ([55]) through Wsii 
given by it is necessary to develop ^{^f cosg, 

(^)3 cos(2A -2/1 + 5) and {^f cos(2A - 2ft. - .g)'^with re- 
spect to the mean anomaly M, the mean longitude of the 
Sun Ls and g, the angle determined in Section [2l Using 
the Kepler's law we obtain this development (See Table[2j 
Table [21 Table 13]). Our value of the eccentricity was taken 
from Simon et al.(1994). 



Table 2. Development of ^(^)'^cos(5) of Venus, t is 
counted in Julian centuries. 



M 


Ls 


g 


period 
d 




cos X 10 " 


tcos xlO ' 








1 


-243.02 




fe^) = 5000344 


-48 


1 





1 


2980.71 




^- ge^) = 50792 


-3582 


1 





-1 


116.75 


(|e- 


+ §e^) = 50792 


-3582 


2 





1 


208.948 




(K) = 516 


-72 


2 





-1 


76.83 




(K) = 516 


-72 


3 





1 


108.27 




%^ = 10 





3 





-1 


,57.25 




f|.^ = 10 






Table 3. Development of ) cos(2(A - ft) - g) of Venus. 
t is counted in Julian centuries. 



M 




g 


period 
d 


cos X 10 ' 


t COS X 10 ^ 





2 


- 1 


76.83 


(1 - |e^) = 9998853 


161 


-1 


2 


-1 


116.75 


(-ie+ jLe-i) = -33859 


2388 


1 


2 


- 1 


57.25 


(|e - !§-e^) = 236993 


-16718 


2 


2 


- 1 


45.62 


(^e^) = 3898 


-550 


-3 


2 


-1 


-2980.71 


48 ^ 





3 


2 


-1 


37.92 


MS 3 ^ 54 


-4 



Table 4. Development of (^) cos(2(A — ft)-|-,g) of Venus. 
t is counted in Julian centuries. 



M 


is 


g 


period 
d 


cos XlO ^ 


(cosxlO ^ 





2 


1 


208.95 


(1 - |e^.) = 9998853 


161 


-1 


2 


1 


2980.71 


(-ie+ -Le^") = -33859 


2388 


1 


2 


1 


108.27 


(|e - i^e-'') = 236993 


-16718 


2 


2 


1 


73.06 


(^e^) = 3898 


550 


-3 


2 


1 


-116.75 


h-' = o 





3 


2 


1 


,55.14 


S#e= = 54 


-4 



The numerical value for Tg is given in Tabled) We 
remind here that our domain of validity is 3000 years as 
in the case in Cottereau and Souchay (2009). 

Thanks to the developments above we give in Table [S] 
and Table [5] the Oppolzer terms, respectively in longitude 
and in obliquity, depending on the dynamical flattening. 
For comparison we give the corresponding nutation coeffi- 
cients (in brackets) of the AMA determined in Cottereau 
and Souchay (2009). They are represented in both FigsH] 
and E] for a 1000 d time span to see the leading oscilla- 
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Table 5. Oppolzer terms in longitude depending on dy- 
namical flattening [A'i's = A/i^: nutation coeflicients of 
the AMA]. 



Argument 


Period 




t sin ujt 


coscjf 




d 


arc second 


irc second/Julian century 


arc second 






(10-7) 


(10-7) 


(10-7) 


2L, 


112.35 


14962988[-21900468] 


-242[-352] 


0[0] 


M 


224.70 


-6134235[889997] 


432651[-62765] 


5416[-61] 


2L. + M 


74.90 


264403 [-346057] 


-18651 [24412] 


-5[-7] 


2L, -M 


224.70 


-76066 [148323] 


5366 [-10461] 


2[-10] 


2Ls + 2M 


56.17 


3466 [-4269] 


-489(602] 


0[0] 


2M 


112.35 


-5749[4521] 


811[-640] 


0[0] 



Table 6. Oppolzer terms in the obliquity depending on 
dynamical flattening [Ae^ — A/^ : nutation coefficients of 
the AMA]. 



rotation g and enter in the denominator during the inte- 
gration are f O'* times larger for the Earth than for Venus. 




200 400 600 800 

time.s(days) 



Fig. 4. The nutation of the AMA (red curve) and the 
nutation of the oppolzer terms (blue and bolt curve) in the 
obliquity of Venus depending on its dynamical flattening 
for a 1000 d span, from J2000.0 



Argument Period cosut t cos tut siiiijj^ 
d arc second arc second/Julian century arc second 
(10~7) (IQ-?) (10"7) 



2L, 


112.35 


-690090 [100741] 


-11[16] 


0[0] 


M 


224.70 


-260831 [0] 


1840[0] 


-248 [0] 


2L, + M 


74.90 


-12184 [15919] 


859[-1123] 


0[0] 


2L, - M 


224.70 


3594[-6822] 


-253 [481] 


0[0] 


2L, + 2M 


56.17 


-160[196] 


22 [-27] 


0[0] 


2U 


112.35 


-122[0] 


17 [0] 


0[0] 



tions. We can remark that the Oppolzer terms are of the 
same order of magnitude as the coefficients of nutation 
themselves. The Oppolzer terms associated with the ar- 
gument M are even larger. This is due to the small value 
of Af+g = -2n/(-243.02/36525)-f 2n/(224.70/36525) = 
76.99 rd/cy which enters in the denominator during the in- 
tegration of the equations of motion. Whereas for the cal- 
culation of the corresponding coefficient of the AMA only 
the numerical value M = 20/(224.70/36525) = 1021.33 
rd/cy appears which is more larger than M + <?. In Figs 
m and O showing the oppolzer terms, the presence of the 
sinusoid with a period of 224 d reflect this fact .Remark 
(table ini) that the Oppolzer terms associated with the ar- 
gument M and 2M have a non zero amplitude, whereas 
the coefficients of nutation of the AMA associated with 
the same arguments do not exist. Indeed we remark here 
that to compute the coefficients of nutation in obliquity 
for this axis, we perform the derivative of W\ with respect 
to h which does not appear in the terms associated with 
the argument M and 2M. Notice also that for the Earth 
(see Kinoshita 1977), the Oppolzer terms depending on 
the dynamical ffattening are negligible with respect to the 
coefficients of nutation of the AMA. The largest Oppolzer 
term in longitude in Kinoshita (1977) is 0" .007559 whereas 
for Venus it is 1" .4962. In obhquity it is 0" .002762 whereas 
for Venus it is 0".6901. The rapid rotation of the Earth 
compared with the slow retrograde rotation of Venus ex- 
plains this contrast : the frequencies which depend on the 




times(days) 



Fig. 5. The nutation of the AMA (red curve) and the 
nutation of the oppolzer terms (blue and bolt curve) in 
longitude of Venus depending on its dynamical flattening 
for a 1000 d span, from J2000.0 



Due to its slow rotation, the triaxiality of Venus 
(1.66 * 10~^) is not negligible compared to the dynamical 
flattening (1.31 * 10~^). Thus the Oppolzer terms depend- 
ing on the triaxiality must be considered and are calcu- 
lated in the following. 

3.3. Oppolzer terms depending on the triaxiality 
Thanks to (gl]), we have : 

Wai ^^^^J[{^fPi{sm6)co,2a]dt 

= -sin^JWaio- ^ psinJ{l + pcosJ)Waii(p) 
^ P=±i 

+ J2 7(l+COSJ)2W^al2(p) (56) 

p=±l 
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where 
W, 



+ isin^/ j {-f cos2{\ - h + l)dt 
+ isin2/ j {'^fcos2{X-h-l)dt 



Waii(p) =i^a'[-^sin2/ j [l^f cos{2pl + g)dt 

— sin /(I — cos/) 
4 

j {-f cos(2A -2h + 2pl + g)dt 

+ - sin /(I + cos/) 
4 

j {'^f cos(2A ~2h- 2pl - g)dt] 
Wai2{p) = Ka'isin^ -^Z^^^)' + 2pg)dt 



(57) 



(58) 



+ -(l-cosI? f (-)3cos2(A- 
4 J r 

+ i(l + cos/)^ 
4 

J {^f cos 2{X - h - pi - g)dt] 



h + pl + g)dt 



and 



3GA/' A-B 



(59) 



(60) 



a3 4 

Using the Hamihon equations and supposing, as for the 
terms depending on the dynamical flattening, that J is a 
small angle, we obtain : 



and : 



-cos/ A^,, + O(J^) 

1 dWal 



2JW^all(l 



(61) 



G sin J dl 
3 . ^dWaw 



1 , ,5I^al 



1 + COS J 



G 



— - COS J ^ 



(1 + C0SJ)2 vdWal2 



4G sin J 



2(P) 



dl 



COS J 



9g 



al2(p) 



(62) 



As is multiplied by (1 — cos J), with our hypoth- 

esis it disappears from the equations (|5T|) and ([5^ . Since 
W^aii(i) and W(j12(p) include / and g in the form oi g + 21 
and 2(7 + 2/ respectively, the last terms in the equation 
((62)) are negligible. Therefore A^j becomes: 



A 



a.J 



2 91^^11(1 
G 5ff 



0(J). 



(63) 



The Oppolzer terms depending on the triaxiality are : 



J sin g , 



2 ,9W^all(l 



sin/ Gsin/ 
+W^aii(i) cos g] + 0(J) 



9g 



■sm5 



A( J cos g) = A J cos g — JAg sin g 

^ dg 



-cosg - Waii(i) sing]. 



(64) 



(65) 



To solve the equations ([M)) and (|55|) it is necessary to 
develop i(^)3 cos(2Z + g), (^)3 cos(2A - 2/1 + 2Z + 5) and 
{-)^ cos{2X—2h—2l—g) with respect to the mean anomaly 
M, the mean longitude of the Sun Lg and the angles I 
and g determined in section (j2.3l) . The coefficients are the 
same as those in Table [3] and S] Only the corresponding 
periods are different, because their calculation includes the 
argument I. The period Ti is very long, as shown in Table 
[TJ Considering our level of accuracy, we suppose in this 
section that e is constant and we take the value of Simon 
et al (1994) as e = 0.006771. 

We can determine the Oppolzer terms depending on 
the triaxiality. 

Table 7. Oppolzer terms in longitude depending on triax- 
iality. Comparison with the corresponding nutation coeffi- 
cients of the AMA in the tables of Cottereau and Souchay 
(2009) 







Oppolzer 


A.M 


Argument 


Period 


sin(cji) 


LC sin{u}t) 




d 


arc second 


arc second 






10-7 


10-7 


2$ 


-121.51 


-11967515 


5994459 


2Ls - 2$ 


58.37 


-3784751 


2880826 


M + 2* 


-264.6 


1490497 


132590 


M -2<S> + 2Ls 


46.34 


-66849 


54201 


M - 2$ 


78.86 


58400 


-39519 


2Ls + 2$ 


1490.35 


5448 


-38866 


-M - 2* + 2Ls 


78.86 


19476 


-13179 


2M + 2<S> 


1490.35 


1062 


-7587 


2M -2<S> + 2Ls 


38.41 


-876 


739 


2M - 2$ 


58.37 


390 


-297 


M + 2<S> + 2Ls 


195.26 


67 


121 


-M + 2<^ + 2Ls 


-264.66 


-263 


-23 


2M + 2<S> + 2Ls 


104.47 


1 


-1 



Table [7] and Table [8] give the Oppolzer terms respec- 
tively in longitude and in obliquity, depending on Venus 
triaxiality. For comparison, we give the corresponding nu- 
tation coefficients of the AMA determined in Cottereau 
and Souchay (2009). They are represented in both Figs 
El and [71 for a 4000 d time span. We remark here also 
that the Oppolzer terms are more important than the 
corresponding coefficients of nutation of the AMA The 
Oppolzer terms associated with the argument 2$ is even 
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Table 8. Oppolzcr terms in obliquity depending on triax- 
iality. Comparison with the corresponding nutation coeffi- 
cients of the AMA in the tables of Cottereau and Souchay 
(2009). 







Oppolzer 


A.M 


Argument 


Period 


cos(a;t) 


CS cos(ajt) 




Q 


arc second 


arc second 






10-7 


10-7 


2$ 


-121.51 


555905 


-275453 


2ls - 2* 


58.37 


-173668 


132365 


M-|-2$ 


-264.6 


-68544 


-6093 


M - 2$ + 2Ls 


46.34 


3074 


2491 


2ls + 2* 


1490.35 


-250 


1786 


M - 2$ 


78.86 


2686 


1816 


-M - 2# + 2Ls 


78.86 


896 


-606 


2M-|-2$ 


1490.35 


-49 


348 


2M - 2$ + 2Ls 


38.41 


-40 


35 


2M - 2$ 


58.37 


18 


14 


M + 2$ + 2Ls 


195.26 


-3 


6 


-M + 2$ + 2Ls 


-264.6 


12 


1 


2M + 2$ + 2Ls 


104.47 









roughly twice larger than the corresponding coefficient of 
the nutation. This is due to the small value oi 21 + g — 
2n/525.81-2n/(243.02/36525) = 944.36 rd/cy which ap- 
pears in the denominator during the integration of the 
equations of motion, whereas for the calculation of the 
corresponding coefficient of the AMA only the frequency 
of the sidereal angle 2$ w 2l~\-2g = 1888.68 rd/cy appears 
which is significantly larger. The appearance of the angle 
g during the integration explains that the other Oppolzer 
terms are more important than the corresponding coeffi- 
cient of the nutation of the AMA, as justified in section 
13.21 for the terms depending on the dynamical flattening. 



1 0-05 



I -0.05 




1000 2000 3000 4000 

times(days) 



Fig. 6. The nutation of the AMA (Blue and bolt curve) 
and the nutation of the Oppolzer terms (red curve) in 
obliquity of Venus depending on its triaxility for a 4000 d 
span, from J2000.0 

Oppositely to Earth, the Oppolzer terms in triaxial- 
ity are not negligible compared to the nominal values of 
the corresponding coefficient of nutation. The frequencies 
depending on the rotation 21 + g which enter in the de- 
nominator during the integration are very small compared 




1000 2000 3000 4000 

tiiiies(days) 



Fig. 7. The nutation of the AMA (Blue and bolt curve) 
and the nutation of the Oppolzer terms (red curve) in 
longitude of Venus depending on its triaxility for a 4000 
d span, from J2000.0 

those of our planet. Now we can determine the nutation 
of the TFA of Venus, which is fundamental from an ob- 
servational point of view. 

4. Numerical results and comparison with the 
motion of the AMA 

Thanks to the equations and we calculate the co- 
efficients of nutation of the TFA. We remind here that the 
nutation is respectively designated in longitude by Ahf 
and in obliquity by AIj . Tables [S] and \TU\ give the coeffi- 
cients of nutation depending on the dynamical flattening 
respectively in longitude and in obliquity. In a similar way. 
Tables [TT] and [T^ give the coefficients depending on the tri- 
axiality. 

Table 9. Avtj^ = A/i/^ : nutation coefficients of the TFA 
in longitude of Venus depending on its dynamical flatten- 
ing 



ArgumenL 


Period 


sin ijL 


t sin ut 


coswt 




d 


arc second 


arc second /Julian century 


arc second 






(10-7) 


(10-7) 


(10-7) 


2L, 


112.35 


-6937480 


-594 





M 


224.70 


-5244238 


369886 


5355 


2L, +M 


74.90 


-81654 


5760 


-2 


2L, - M 


224.70 


72257 


-5095 


-8 


2L, + 2M 


56.17 


-803 


113 


0[0] 


2M 


112.35 


-1228 


171 


0[0] 



In this section, we will show the difference between 
the nutation of the TFA, calculated in this paper and 
that of the AMA of Venus {51 and A/i), as calculated 
by Cottereau and Souchay (2009). Figs [8] and [9] represent 
respectively the nutation in longitude and in obliquity of 
the two axes for a 4000 d time span. 

In longitude we can point out two important speciflc 
remarks : 

— The nutation of the TFA is significantly smaller than 
the nutation of the AMA. Indeed the amplitude peak 
to peak, of the nutation of the TFA (an amplitude 
of 1.5") is twice smaller than that of the AMA (an 
amplitude of 3" ) 
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Table 10. Ae/s — ^^fs '■ nutation coefficients of the TFA 
in obliquity of Venus depending on its dynamical flatten- 
ing 



Argument 


Period 


cosa;i 


I COSLjf 


sin u?f 




d 


arc second 


arc second/Julian century 


arc second 






(10-7) 


(10-7) 


(10-7) 


2La 


112.35 


-589348 


5 





M 


224.70 


-260831 


1839 


-248 


IL, + M 


74.90 


3735 


-264 





1L, - M 


224.70 


-3228 


-228 





2L. + 2A/ 


56.17 


37 


-4 





2M 


112.35 


-122 


17 






Table 11. fa = A/i/a nutation coefficients of the TFA 
in longitude of Venus depending on its triaxialit 



Argument 


Period 


sin(:(jf) 




d 


arc second 
10-7 


2<I> 


-121.51 


-5973056 


ILs - 2$ 


58.37 


-903925 


Af + 2* 


-264.6 


1623087 


M - 2* -1 2Ls 


46.34 


-12648 


i\/-2S 


78.86 


29199 


2Ls + 24> 


1490.35 


-33418 


-M - 25 + 2Ls 


78.86 


6297 


2M - 2# -1 2Ls 


38.41 


-137] 


2M - 2# 


58.37 


195 


M + 2<S> + 2Lb 


195.26 


188 


-M + 2S + 2Ls 


-264.66 


-286 


2M + 2# -1 2Ls 


104.47 







1000 2000 3000 4000 



Fig. 8. The nutation of the TFA (blue curve) and the 
nutation of the momentum axis in longitude of Venus for 
4000 d time span, from J2000.0 




1000 2000 3000 4000 

times(days) 



Fig. 9. The nutation of the TFA (blue curve) and the 
nutation of the momentum axis in obliquity of Venus for 
4000 d time span, from J2000.0 



Table 12. Ae/a = A//a : nutation coefficients of the TFA 
in obliquity of Venus depending on its triaxiality 



Argument Period cos(ujt) 

d arc second 



2<I> -121.51 280452 

2ls - 2$ 58.37 -41303 

M + 2S -264.6 -74637 

M -2<i> + 2Ls 46.34 5565 

2ls + 2S 1490.35 1536 

M - 2* 78.86 4502 

-M - 25 + 2Ls 78.86 290 

2M + 2# 1490.35 299 

2M - 2# -1 2Ls 38.41 -5 

2M - 2* 58.37 32 

M + 29> + 2Ls 195.26 3 

-M + 2$ + 2Ls -264.6 13 

2M + 2l> + 2Ls 104.47 



— The nutation of the TFA is dominated by three sinu- 
soids associated with the arguments 2Ls, M and 2$, 
with respective periods 112.35 d, 224.70 d and 121.51 
d whereas the nutation of the AMA is dominated by 
two sinusoids of argument 2Ls and 2$. 

The same kind of remark is available in obliquity but the 
difference between the nutations is relatively less impor- 
tant. The amplitude of the nutation of the TFA is varying 
between 0".10 and 0".08 (peak to peak) whereas the am- 
plitude of the AMA is varying between -Q."14 and 0."11. 
The nutation in obliquity is also dominated by a sinusiod, 



with a period M which is not as large as the corresponding 
sinusoid of the AMA. 

Finally we can highlight the differences between the 
Earth and Venus. For the Earth, the nutation of the 
two axes (angular momentum and third figure axis) 
are roughly the same (Woolard,1953, Kinoshita, 1977) 
whereas in the case of Venus they are significantly dif- 
ferent. We can also remark that the leading component of 
nutation of the third Venus figure axis in longitude due 
to the gravitational action of the Sun, with argument 2Ls 
(seeO has an amplitude of 0" .693. This is of the same or- 
der as the leading 2Ls the nutation amplitude of the Earth 
due to the Sun, i.e. 0".998 despite the fact that Venus 
has a very small dynamical flattening. As explained by 
Cottereau and Souchay (2009) this is due to the compen- 
sating role of the very slow rotation of Venus. Moreover, 
notice that in the case of Venus the argument Ls stands 
for the longitude of the Sun as seen from the planet, so 
that the corresponding period of the leading nutation term 
with 2Ls argument is 112.35 d, whereas it is 182.5 d in the 
case of the Earth. 

5. Determination of the indirect planetary effects 
on the nutation of Venus 

Using the ephemeris DE405 (Standish, 1998), we have 
computed the nutation of Venus by numerical integration 
with a Runge-Kutta 12*'* order algorithm. Figs [TU] and 
111! show the differences between the nutations in obliq- 
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uity and in longitude of the AMA as computed from the 
analytical tables (Cottereau and Souchay, 2009) and that 
from the numerical integration for a 4000 d time span. 
The residuals obtained consist clearly in periodic com- 
ponents with small amplitudes, of the order of 10~^" in 
obliquity and 10"'^" in longitude. This numerical integra- 
tion validates the results of Cottereau and Souchay (2009), 
down to a relative accuracy of 10~^. Moreover Kinoshita's 
model used in Cottereau and Souchay (2009) supposed 
a Keplerian motion of Venus around the Sun. It is well 
known that the effects of planetary attraction into Earth's 
orbit (called indirect planetary effect) entails a departure 
from the Keplerian motion and that this departure in- 
duces new nutation terms as calculated by Souchay and 
Kinoshita (1996). In order to infer whether the discrepan- 
cies between our numerical integration and our analytical 
computation are caused by this indirect planetary effect 
in Venus orbit, we performed a spectral analysis of the 
residuals. 

0.00005 ^ — . — , — . — , — . — . — . — . — , — , — , — . — . — , — , — . — . — ^ 




-0.00002 



1000 2000 3000 4000 

temps(days) 

Fig. 10. Difference between the nutation in obliquity of 
the AMA and the numerical integration for 4000 d time 
span, from J2000.0. The curve at bottom represents the 
residual after substracting the sinusoidal terms of the 
Table [H 



0.0020 - 




1000 2000 3000 4000 

temps(days) 



Fig. 11. Difference between the nutation in longitude of 
the AMA and the numerical integration for 4000 d time 
span, from J2000.0. The curve at bottom represents the 
residual after substracting the sinusoidal terms of the 
Table m 

The leading oscillations of the two signals (in longitude 
and in obliquity) are determined thanks to a fast Fourier 
Transform (FFT). Tables[T51and[Hlgive the leading ampli- 
tudes and periods of the sinusoids characterizing the signal 



Table 13. Coefhcient of rigid Venus nutation from the 
indirect planetary contribution in longitude. Comparison 
with the respective value in tables of Souchay and 
Kinoshita (1995) when we have the Earth as a indirect 
planetary contribution. 



Planetary effects 


Period 


amplitude 


Amplitude 


SK 


SK 






sin 


cos 


sin 


cos 




.yr 


O.Olmas 


O.Olmas 


O.Olmas 


O.Olmas 


2V-2E 


0.80 


31.5 


-9.6 


-9.6 


0.0 


2V-3M 


0.60 


19.3 


4.0 


/ 


/ 


3V-3E 


0.53 


-2.0 


13.4 


0.5 


-0.2 


V-,J 


0.64 


15.1 


-26.2 


/ 


/ 


V-E 


1.60 


-2.6 


-15.3 


6.6 


0.0 




4.11 


-13.3 


-38.4 


/ 


/ 


2V+2M+2,J 


0.22 


6.2 


12.7 


/ 


/ 




9.31 


-274.5 


421.0 


/ 


/ 




0.20 


-1.5 


-2.3 


/ 


/ 




0.70 


-21.0 


-6.5 


/ 


/ 


2V-2Me 


0.19 


8.5 


-7.1 


/ 


/ 


2V+J+M 


0.25 


4.6 


-2.9 


/ 


/ 



Table 14. Coefhcient of rigid Venus nutation from the 
indirect planetary contribution in obliquity. Comparison 
with the respective value in tables of Souchay and 
Kinoshita (1995) when we have the Earth as a indirect 
planetary contribution. 



Planetary effects 


Period 


amplitude 


amplitude 


SK 


SK 






sin 


cos 


sin 


cos 




yr 


O.Olmas 


O.Olmas 


O.Olmas 


O.Olmas 


2V4-2M-I-2J 


0.22 


0.5 


-0.4 


/ 


/ 


2S 


14.96 


0.7 


0.4 


/ 


/ 




4.06 


-0.3 


0.2 


/ 


/ 


2V-2Mer 


0.19 


-0.3 


-0.3 


/ 


/ 


V-2.J+M 


0.50 


-0.3 


0.0 


/ 


/ 



in FigsfTUlandfTTjwhere the curve at the bottom represents 
the residuals after subtraction of these sinusoids. The pe- 
riods presented in the Tables [T^ and ITU do not correspond 
to any period of the tables given in the precedent section 
starting from the keplerian approximation. On the oppo- 
site, when comparing them with the tables of nutation of 
the Earth taken from Souchay and Kinoshita (1996), simi- 
lar periods appear which correspond to the combination of 
planetary longitudes. This is a clear confirmation that the 
differences between our analytical computation and our 
numerical integration, in Figs [TU] and 1111 are essentially 
due to the indirect planetary effects, negligible at first or- 
der and not taken previously into account by Cottereau 
and Souchay (2009). Tables flSl and [T4l also present, when 
they have been clearly identified, the combination of the 
planetary longitudes corresponding to the detected sinu- 
soids. Moreover, we give, when they are available, the cor- 
responding amplitudes of the nutation of the Earth due 
to the indirect effect of Venus. We can thus point out the 
similitude and reciprocity of the indirect planetary effects 
of Venus on the rotation of the Earth, and of the indirect 
planetary effects of the Earth on the rotation of Venus. 
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6. Conclusion and Prospects 

In this paper we achieved the accurate study of the ro- 
tation of Venus, for a rigid model and on a short time 
scale begun by Cottereau and Souchay (2009), by apply- 
ing analytical formalisms already used for the rigid Earth 
(Kinoshita, 1972, 1977). The differences between the rota- 
tional characteristics of Venus and our planet, due to the 
slow rotation of Venus and its small obliquity, have been 
highlighted. 

Firstly we have precisely determined the polhody, i.e 
the torque free rotational motion for a rigid Venus. We 
have adopted the theory used by Kinoshita (1972, 1992) 
and we have given the parametrization and the equations 
of motion to solve the motion. We have shown that the pol- 
hody is significantly elliptic, on the contrary of the Earth 
for which it can be considered as circular in first approx- 
imation. Moreover it is considerably slower. Indeed, the 
period of the torque free motion is 525.81 cy for Venus 
whereas it is 303 d for our planet, when considered as 
rigid. 

Then we have determined the motion of the third fig- 
ure axis, which is fundamental from an observational point 
of view. We have calculated the Oppolzer terms due to the 
gravitational action of the Sun using the equation of mo- 
tion of Kinoshita (1977) as well as the corresponding de- 
velopment of the disturbing functions. We have compared 
them with the coefficients of nutation for the angular 
momentum axis taken in Cottereau and Souchay (2009). 
One of the important results is that these Oppolzer terms 
depending on the dynamical flattening are of the same 
order of amplitude as the coefficients of nutation them- 
selves, whereas for the Earth (Woolard, 1953; Kinoshita, 
1977) these Oppolzer terms are very small with respect 
to the coefficients of nutation of the angular momentum 
axis. Moreover we have computed the Oppolzer terms de- 
pending on the triaxiality, which is not done in Kinoshita 
(1977) in the case of the Earth for which they are negligi- 
ble. In the case of Venus, these Oppolzer terms are signif- 
icant. Even larger than the corresponding coefficients of 
nutation of the angular momentum axis. 

Thanks to our Oppolzer terms we have also been able 
to give the tables of nutation of the third flgure axis from 
which we computed the nutation for a 4000 d time span. 
The comparison with the nutation of the angular momen- 
tum axis, calculated from Cottereau and Souchay (2009), 
is also given. The nutation of the third figure axis is sig- 
nificantly smaller peak to peak than the nutation of the 
angular momentum axis in longitude, and less important 
in obliquity. The amplitude of the largest nutation coef- 
ficient in longitude of the third figure axis (1.5") is half 
the one of the angular momentum axis (3"). The ampli- 
tude of the nutation in obliquity of the third figure axis is 
0".18 peak to peak whereas the amplitude of the angular 
momentum axis is 0" .25. The nutations of the third fig- 
ure axis, in obliquity and in longitiide, are dominated by 
three sinusoids associated with the arguments 2Ls, M and 
2$, with respective periods 112.35 d, 224.70 d and 121.51 



d. The nutation of the angular momentum axis is dom- 
inated by two sinusoids with the argument 2Ls and 2$. 
Our results have shown that although the axis of angu- 
lar momentum and the third figure axis can be considered 
identical in the case of the Earth (Kinoshita, 1977), this 
approximation does not hold in the case of a slowly rotat- 
ing planet as Venus. 

At last we have validated our analytical results down to 
a relative accuracy of 10^^ with a numerical integration. 
Moreover, we have confirmed, by using results in Souchay 
and Kinoshita (1996) for the nutation of a rigid Earth, 
that the differences between our analytical compiitation 
and our numerical integration are essentially due to the 
indirect planetary effects, which was not taken into ac- 
count by Cottereau and Souchay (2009). We think that 
this study is fundamental to understand the behavior of 
Venus rotation in a very accurate and exhaustive way for 
short time scales, and should be a necessary starting point 
to another similar study including non rigid effects (elas- 
ticity, atmospheric forcing etc.). 
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7. appendix 

7.1. Development for a small value of the triaxiality 
(Kinoshita, 1972) 

We note b - 



I = [- l)esin2[ 

+{^(b'' + 6b+l)e\md+0{e^), 



(66) 



9 =9 + 



1 



-esin2Z 



--^(6+ 1)6^ sind 



(67) 



J ^ J ^ (26+ 1)6^ tan J 

16 

+ tan J(^ecos2[- ^6e^ cos4Z~) + O(e^), 
with : 



(68) 



~ 1 13 

J = j+2^ tan j + tan j(- + — tan^ j) + 0{e^) 

G 1 

l = nixt with ni = — cos J[l - - (6^ + 3)e^] + O(e^) 
U 8 

g = rig X t with 



7.2. Development for a small value of the angle J 
(Kinoshita, 1972) 

The polar angles I and J leading to the determination of 
the free rotational motion are given by : 



4 y 1-e' l + ecos2Z ' 



+1(1 + e)f [-^^^ + ^=l^(r - 0] ) + 0(/) 

4 l + ecos2t vl — ^ 

J = jyi + Y^cos2[+0(j3) (70) 



with 



tanT =\- — ^tanf 
V 1 + e 



'■g 



'(T + 4)G+^(i-^r^) 



2M D' 



x[l + y[^i2] + 0(/). (71) 



